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Abstract 

The /3-decay half-lives of neutron-rich nuclei with 20 ^ Z ^ 50 are systematically investigated 
using the newly developed fully self-consistent proton-neutron quasiparticle random phase approxi- 
mation (QRPA), based on the relativistic Hartree-Fock-Bogoliubov (RHFB) framework. Available 
data are reproduced by including an isospin-dependent proton-neutron pairing interaction in the 
isoscalar channel of the RHFB-I-QRPA model. With the calculated /3-decay half-lives of neutron- 
rich nuclei, a remarkable speeding up of r-matter flow is predicted. This leads to enhanced r-process 
abundances of elements with A > 140, an important result for the understanding of the origin of 
heavy elements in the universe. 
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Nuclear /3-decay plays an important role not only in nuclear physics but also in other 
branches of science, notably astrophysics and particle physics. In nuclei the investigation 
of /3-decay provides information on the spin and isospin dependence of the effective nuclear 
interaction, as well as on nuclear properties such as masses Q, shapes 3, and energy 
levels |3|] . In nuclear astrophysics /3-decay half-lives set the time scale of the rapid neutron- 
capture process (r-process), and hence determine the production of heavy elements in the 
universe j4-6|. In particle physics /3-decay was used to obtain the first experimental evidence 
of parity violation and can be utilized to verify the unitarity of Cabibbo-Kobayashi- 
Maskawa (CKM) matrix 9|. 

Important advances in the measurement of nuclear /3-decay half-lives have been achieved 
in recent years with the development of radioactive ion-beam facilities, especially for nuclei 
near the neutron shell closures of 50 and 82 10Nl2|. Quite recently /3-decay half-lives of 



very neutron-rich Kr to Tc isotopes with neutron number between = 50 and 82 have been 
measured j3] • The new experimental results indicate a systematic deviation from half-lives 
predicted by standard calculations based on the finite-range droplet model (FRDM) plus 

y measured (3- 



14j . and it has 



quasiparticle random phase approximation (QRPA). The effects of the new 
decay half-lives on r-process nucleosynthesis have been investigated in Ref. 
been shown that the main effect is an enhancement in the abundances of isotopes with mass 
number A = 110 — 120, relative to abundances calculated using /3-decay half-lives estimated 
with the FRDM+QRPA. 

Most neutron-rich nuclei relevant for the r-process are still out of the reach of experimental 
facilities and, therefore, r-process calculations are based on theoretical predictions for /3- 
decay half-lives. Theoretical investigations of the nuclear /3-decay started in the 1930's, 



with the introduction of the famous Fermi theory [15|]. Two types of microscopic approaches 
are nowadays mainly used in large-scale calculations of nuclear /3-decay half-lives, that is , the 
shell model and the proton-neutron QRPA. Specifically the shell model has been applied 
to calculation of /3-decav half-lives for nuclei at A^ = 50, 82, 126 and reproduces well the 
experimental half-lives {g, | ' 
from the magic numbers are not feasible yet because of extremely large configuration spaces. 
Compared to this approach, the proton-neutron QRPA can be applied to arbitrary heavy 
systems. Nuclear /3-decay calculations have been carried out using the QRPA based on the 
FRDM [isl, the extended Thomas- Fermi plus Strutinsky integral (ETFSI) model 19|, the 



f-lives for nuclei at A^ 

However, shell-model calculations for heavy nuclei away 
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Skyrme Hartree-Fock-Bogoliubov (SHFB) model [20|], and the density functional of Fayans 
(DF) |21|. 

Covariant density functional theory (CDFT) has been app lied very successfully to the 



22-I2J]. The CDFT framework 



description of a variety of nuclear structure phenomena 
naturally includes the nucleon spin degree of freedom, and the resulting nuclear spin-orbit 
potential automatically emerges with the empirical strength, thus producing a good agree- 
ment with the experimental single-nucleon spectrum 25||. Based on the relativistic Hartree- 
Bogoliubov (RHB) model in the CDFT framework, the QRPA has been formulated 26| and 
employed in the calculations of /9-decay half-lives of neutron-rich nuclei in the regions of 
N ^50 and N ^ 82 



27 



28 



To reliably predict properties of thousands of unknown nuclei relevant to the r-process, 
the self-consistency of the QRPA approach is essential. Only recently the fully self-consistent 
relativistic RPA has been formulated based on the relativistic Hartree-Fock (RHF) the- 
ory 29[. Using the RHF+RPA model excellent agreement with data on the Gamow- Teller 
(GTR) and spin-dipole resonances in doubly magic nuclei has been obtained, without any 
readjustment of the parameters of the covariant energy density functional 29|, |30|. Recently 
also the relativistic Hartree-Fock-Bogoliubov (RHFB) theory has been developed, providing 



31 



32|. 



a unified description of both mean field and pairing correlations 

In this Letter we report the implementation of a fully self-consistent proton-neutron 
QRPA based on the RHFB theory, and its first systematic application in calculations of /3- 
decay half-lives of neutron-rich nuclei with 20 ^ Z ^ 50, extending over the whole r-process 
path from = 50 to = 82. The effect on r-process nucleosynthesis simulations is also 
investigated using the classical r-process model. 



The details of the QRPA formalism in the canonical basis can be found in Refs. [20|, |26 |. 
In the RHFB+QRPA model both the direct and exchange terms are taken into account, so 
the matrix elements of the particle-hole (p-h) V^^ and particle-particle {p-p) V^^ residual 
interactions are respectively written as 



C'n' = / j dr,dr,f;{r,)f^,{r,)Y,V^X^,2) 

[/p'(^2)/n(n) - fn{r2)fp'{ri)] , 



(1) 
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J J T=1,0 
[/n'(^2)/p'(n) - /p'(^2)/n'(n)] • 



(2) 



p,p', and n,n' denote proton and neutron quasiparticle canonical states, respectively, and (pi 
and T denote corresponding coupling channels. /p(„) are canonical wave functions extracted 
from the RHFB calculations. 

Because of the presence of exchange terms, the proton-neutron QRPA interaction contains 
terms generated not only by the isovector meson exchange (p and vr), but also by the 
exchange of isoscalar mesons {a and u). As in Ref. 29|] the pionic zero-range counter term 
introduced to remove the contact part of the pseudovector vr-N coupling reads 



-T7o757 



-T7o757 



5(n - 



(3) 



where the self-consistency of the model requires o' = 1/3. For the isovector (T = 1) p-p 
channel the pairing part of the Gogny force DIS 33| is consistently used both in the RHFB 
and QRPA calculations. For the isoscalar (T = 0) proton-neutron pairing in the QRPA we 



employ a similar interaction that was previously used in Refs. 

2 



20 



28| 



^T=o(l,2) 



-[{rl-r2)/^lj] 



n 



5=1,T=0 



(4) 



with /ii = 1.2 fm, /i2 = 0.7 fm, gi = 1, g2 = —2. The operator ns=i t=o Projects onto states 
with 5 = 1 and T = 0. Vq is the overall strength of the T = proton-neutron pairing. 

With the individual transition strengths Bm obtained from QRPA calculations, the (3- 
decay half-life of an even-even nucleus is calculated in the allowed Gamow- Teller approxi- 
mation using the expression: 

D 



(5) 



where D = 6163.4 ± 3.8 s and gA = 1- The sum runs over all final states with an excitation 
energy smaller than the value. The integrated phase volume f{Z, Em) can be written as 

"Em 



f{Z,E„ 



PeEe{Em — Eg) Fo{Z , Ee)dEe, 



(6) 



where Pe, E^, and Fq{Z,E(,) denote the emitted electron momentum, energy, and Fermi 



function, respectively 



6|. The /9-decay transition energy E^, that is, the energy difference 
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FIG. 1: (Color online) Nuclear /3-decay half-lives of Fe (left) and Cd (right) isotopes, calculated 



10|]. Open tri- 



with the PKOl effective interaction [34!], compared with the experimental values 
angles, diamonds, and squares denote values obtained using the RHFB+QRPA with the strength 
parameter of the T = pairing: Vq = 0, 155, and 255 MeV, respectively. The RHFB+QRPA values 
obtained with the Vq of Eq. dH) are denoted by open circles. 

between the initial and final state, is calculated using the QRPA: 

Em. = ^np — -E'qRPA, (7) 

where -Eqrpa is the QRPA energy with respect to the ground-state of the parent nucleus and 
corrected by the difference of the neutron and proton Fermi energies in the parent nucleus 



20l |. and A„p is the neutron-proton mass difference. 

Nuclear /3-decay half-lives are very sensitive to the T = proton-neutron pairing in- 
teraction, and its strength Vq is determined by adjusting QRPA results to empirical half- 



lives 



20 



27 



28|. Taking '''^Fe and ^^'^Cd as reference nuclei for the two mass regions, the 
value of Vq is determined as 155 MeV and 255 MeV, respectively. Using these two values, 
the calculated half-lives of Fe and Cd isotopes are shown in Fig. [1] For comparison, the 
experimental values and the results of a calculation without the T = pairing are also 
displayed. One notices that the /3-decay half-lives calculated without the inclusion of T = 
pairing are systematically much longer than the experimental half-lives, both for Fe and Cd 
isotopes. In previous studies a constant value of Vq was taken for one isotopic chain or one 
mass region, determined by adjusting to the known half-lives of selected nuclei in the isotopic 
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FIG. 2: (Color online) Contour map of /3-decay half-lives for the Z = 20 
The experimental half-lives Q, llj ] and the RHFB-I-QRPA results obtained with the effective 
interaction PKOl are shown in panel (a) and panel (b), respectively. For reference the r-process 
path is also displayed. 



chain [27|, or several nuclei in the same mass region [20|. This procedure, of course, limits 
the predictability of the model. Moreover, as shown in Fig. [1] for the Cd isotopes, when 
Vq is determined by the /3-decay half-life of ^"^^Cd the calculated results underestimate the 
half-lives of ii8'12o,i22q^ This indicates that the half-lives of an isotopic chain cannot always 
be reproduced using a constant value Vq, and points out to a possible isospin-dependence of 
Vq. In the present work, therefore, we employ the following ansatz for an isospin-dependent 
pairing strength: 

Vd 



Vq = Vl + 



(8) 



1 _^ ga+6(7V-Z) ' 

and adjust the parameters to reproduce the known half-lives of even-even nuclei with 20 ^ 
Z ^ 50: Vl = 134.0 MeV, Vd = 121.1 MeV, a = 8.5, and b = -0.4. 

With the isospin-dependent strength Eq. (|8]) of the proton-neutron pairing interaction, 
the calculated /3-decay half-lives of both the Fe and Cd isotopic chains are in excellent 
agreement with data. Encouraged by this success, we proceed to calculate the half-lives 
of even-even nuclei with 20 ^ Z ^ 50 using the RHFB+QRPA model, and compare the 
theoretical values to data in Fig. [2l Even though a wealth of new data on nuclear /3- 
decay half-lives have been obtained in recent years, only few measurements can reach the 
r-process path, especially for the r-process path around = 82. The present RHFB+QRPA 
calculation yields results in good agreement with the data, especially for nuclei with half-lives 
Ti/2 < 1 s. Only an overestimation of half-lives for Ni, Zn, and Ge isotopes near stability 
line is noticed. The longer half-lives predicted for Ni isotopes is a common problem in self- 



consistent relativistic QRPA calculations 



28|. For Zn and Ge isotopes, the differences 




FIG. 3: (Color online) Nuclear /3-decay half-lives of Cd isotopes calculated in the RHFB+QRPA 
framework using the isospin-dependent proton-neutron pairing interaction, and the effective in- 
teraction PKOl. For comparison, the experimental values 



10l |. as well as the theoretical results 
0,the SHFB+QRPA 



20(, and the FRDM+QRPA calculation 



181], are 



obtained in the RHB+QRPA model with the effective interaction DD-MEl 
calculation with the effective interaction SkO' 
also shown. 

between present results and the experimental values are remarkably reduced as the neutron 
number increases. Taking ®^Ge for example, the theoretical result is 1.3 s, rather close to 
the experimental value: 0.954 ± 0.014 s. 

In Fig. 131 the calculated /3-decay half-lives of a sequence of Cd isotopes using 
RHFB+QRPA approach are displayed in comparison with measured values and previous 
theoretical results. One notices that the RHFB+QRPA model reproduces in detail the em- 



pirical /3-decay half-lives, whereas the RHB+QRPA and SHFB+QRPA (20| calculations 
that used a constant proton-neutron pairing strength cannot yield the appropriate isospin 



dependence of /3-decay half-lives. The global FRDM+QRPA calculation [18| systematically 
overestimates the measured half-lives of Cd isotopes. It has been pointed out that the over- 
estimation of half-lives in the FRDM+QRPA approach can be at least partially attributed 



to the omission of the T = pairing 20[. This is further confirmed in this work as the 
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FIG. 4: (Color online) The impact of nuclear /3-decay half-lives on the calculated r-process 
abundance. The solid (dashed) curves correspond to r-process abundance calculated with the 
RHFB-hQRPA (FRDM-FQRPA) /3-decay half-lives in comparison with the data denoted by the 



points. In all calculations, nuclear masses are taken from the mass evaluation 



lOl if available, 



36l | are employed. Panels (a)-(d) correspond to the 



otherwise predictions of RMF mass model 
times Tr = 1.5,2.0,2.5, and 3.0 s, respectively. 

half-lives are systematically reduced with the inclusion of T = pairing. Furthermore, 
for the Cd isotopes with ^ 90 the FRDM+QRPA calculation predicts shorter half-lives 
compared to the RHFB+QRPA. In the FRDM a shape transition occurs from the spherical 
^^^Cd to the deformed ^^^Cd with the quadrupole deformation P2 = 0.125. Therefore, the 
shorter half-lives for Cd isotopes with ^ 90 in FRDM+QRPA calculation may be due 
to the effect of nuclear deformation. Similar systematics are also found for the Zr, Mo, 
Ru, Pd isotopes where shape transitions occur. However, it should be noted that a recent 
self-consistent QRPA calculation that used a Skyrme interaction has shown the opposite 



trend [35 1 



To analyze the impact of the predicted /3-decay half-lives on r-process abundances, we 



37 



with 



have also performed a classical r-process calculation similar to those in Refs. 
neutron density n„ = 10^^ — 10^^ cm~'^ and temperature T = 1.5 x 10^ K. Figure S] displays 
snapshots of r-process abundance for different process time r. It is shown that the half-lives 
calculated with RHFB+QRPA model (in solid lines) produce a faster r-matter flow at the 
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= 82 region, thus yield higher r-process abundances of elements with A > 140. By 
summing up the half-lives of r-path nuclei at = 82, one can roughly estimate the time 
when the r-process passes the = 82 shell. Based on the RHFB+QRPA results, this time 
is speeded up to 0.25 s from the 1.29 s predicted by the FRDM-I-QRPA calculation. This is 
an important result for the estimate of the duration of the r-process, and hence the origin 
of heavy elements in the universe. In addition, it is found that the abundance at A ~ 80 
is higher using the RHFB+QRPA result, a result that can be easily understood because it 
takes more time to pass the r-path nuclei ^^Ni and ^°Zn and, as a result, higher abundances 
are accumulated. 

In summary, this work introduces the fully self-consistent proton-neutron quasiparticle 
random phase approximation (QRPA), based on the relativistic Hartree-Fock-Bogoliubov 
(RHFB) framework. Using an isospin-dependent proton-neutron T = pairing interaction, 
the RHFB+QRPA model is applied to a calculation of /3-decay half-lives of neutron-rich 
nuclei with 20 ^ Z ^ 50, extending over the entire r-process path from A^ = 50 to A^ = 82. 
It is found that the RHFB+QRPA model calculation reproduces the experimental /3-decay 
half-lives for neutron-rich nuclei, especially for nuclei with half-lives less than one second. 
Using the calculated /3-decay half-lives of neutron-rich nuclei, a remarkable speeding up of r- 
matter flow is predicted. This leads to an enhancement of r-process abundances of elements 
with A > 140. 
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